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Abstract 

In this paper we address again the problem of the connection between mul- 
titime Brownian sheet and heat type PDEs. The main results include: the volu- 
metric character of the solutions of the forward (backward) diffusion-like PDEs; 
the forward mean value of a Brownian process as the solution of the forward heat 
PDE; the backward mean value of a Brownian process as the solution of the back- 
ward heat PDE; the multitime stochastic processes with volumetric dependence; 
the stochastic partial derivative of a stochastic process with respect to a multitime 
Wiener process; Hermite polynomials stochastic processes; union of Tzitzeica hy- 
persurfaces (constant level sets of multitime stochastic processes with volumetric 
dependence). The original results permit to extend the complete integrability the- 
ory to multitime stochastic differential systems, using path independent curvilinear 
integrals and volumetric dependence. 



1 Introduction 

The papers lIS have studied stochastic curvilinear integrals along all sufficiently 
smooth curves in R^. The most simple situation is those of increasing curves. 

Our papers ||9l, US, 1 14] have extended this point of view to stochastic curvilinear 
integrals in R™ and to completely integrable stochastic differential systems. These 
research trends and the original results are based on It6-Udri§te stochastic calculus 
rules 

dWf dW'^ = 5"'' Ca(t)df, dW} dt" = dt" dW^ = 0, df dt^ = 0, 

where a,b - l,d; a,/5 = l,m, the parameter f = (f', f"') E W" means the multitime, 
6"'' is the Kronecker symbol, v — f ^ • ■ • f'" is the volume of the hyperparallelepiped 
Qor c R"_l, Ca(t) - ^ and the tensorial product d"'^ c„{t) represents the correlation 
coefficients. These formulas can be extended to stochastic processes on Riemannian 
manifolds or on manifolds with a fundamental tensor of type (2, 0) that determines the 
correlation coefficients. A similar single-time theory can be found in the paper |[T2l . 
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Also some ideas have the roots in the basic book f4l containing the single-time theory 
of Brownian motion, stochastic integral and differential stochastic equations. 

The main objective of the paper is the connection between multitime Brownian 
sheet and heat type PDEs. Section 2 shows that the solutions of the forward diffusion- 
like PDEs system have a volumetric character and that the forward mean value of a 
Brownian process is the solution of the forward heat PDE. Section 3 proves that the 
solutions of the backward diffusion-like PDEs system have a volumetric character and 
that the backward mean value of a Brownian process is the solution of the backward 
heat PDE. Section 4 describes the differentiable multitime stochastic processes as sums 
of series of Hermite polynomials. Section 5 studies unions of Tzitzeica hypersurfaces 
and their connection to differentiable multitime stochastic processes. Some related 
deterministic ideas can be found in the papers Q, ifTOl . ifTTl . 



2 Forward diffusion-like PDEs system 

Let us now consider a real-valued function / (f, x) , t - (t^, f") e R'", x - (x^, x^) e 
R'', with f(t, 0) — 0, which has continuous partial derivatives of the first order with re- 
spect to t" , and of the second order in x" . If 

W, = (W,"), t = (f', f'") 6 Rf, fl = M 

is a multitime vectorial Wiener process |9|, then the function / defines a stochastic 
process f{t,Wt), t e R™. Let us accept (see |1|) that /(f, W,) is a solution of the 
forward diffusion-like PDEs system 

^ it, W,) - i cjt) AJ(t, W,) = 0, for a^T^, (2. 1) 

where 

dx"dx" 

Theorem 2.1 The solution f of the forward diffusion-like PDEs system (1) depends on 
the point (f', f") only through the product of components v — f' • ■ • f", i.e., it is a 
function of the volume v — t^ ■ ■ ■ f" o/Qor C R™. 

Proof. From (1) it follows 

c/i(t) ^ it, W,) = c„(f) ^ (f, W,),a+ p. 

The general solution of this PDEs system is f{t) - ip(t^ ■ ■ ■ f\ W,), t e R"'. 

If V = f' • • ■ f'" is the volume of Qq, c R'" and g (v, x) ip[t^ ■ ■ ■ f™, has continu- 
ous partial derivatives of the first order in v and of the second order in x, then from (1) 
it follows that g satisfies t\\t forward multitime heat PDE 

^ - ^Ag = 0. (2.2) 
ov 2 
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Lemma 2.2 (see lIU) Let 



W:^ ^{Wy, t = it\...,f")eR'^,xeR\a= l,d 

be a family of Brownian sheets with P{W^ — x] — \ for t e dR'". If f : R'' ^ R is a 
continuous bounded function and u(t, x) — E[f(W^)] is the forward mean value, then 
u(t, x) is a unique bounded solution of forward heat type PDEs 



du 1 
— (t,x) - - 
dt" 2 



(f,x) = -Ca{t)l^Mt,x), t eR'"; 



u(t,x)^fix), tedR". 

Proof. . Let K(t; x, y) be the transition density of an m-parameter valued Brownian 
sheet Wit) - W^(f', ...,?'"), starting on dR"^ from x € R'', under P, and moving to 
y € R'' at muhitime f = it\..., f) e R"^, i.e., 

1 1 /-llx-vlH 
K(t;x,y)^—P[W-^(t)edy] = —- — 77. exp " ^" 



dy (Int^ ■■■t"'yi^ ^\2f'---/ 

This function is a solution of the forward heat type PDEs, i.e., 

d 1 

—K{t;x,y) = -Ca^Mnx,y). 

Taking into account the definition of the mean value u{ t, x), the assumptions on the 
function / and the differentiability of the improper integral function, we find 



£ = i = lmlmx,y)dy 



\ca \ f(y)A,,K(t;x,y)dy 

N't 
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1 r 1 

-Ca A., f(y) K(t; x, y) dy ^ -c^ A^u. 

2 Jr'I 



Jr-i 2 
By definition and by the continuity and boundedness of the function / we easily have 

lim u(s,y) - u{t,x) - fix), 

A' 

for f e dR"^ and x e R''. The existence and uniqueness Theorem for a Cauchy problem 
associated to heat PDE ends the proof. ■ 

The mean value function u(t, x) is a function of the volume v = f' • ■ ■t"\ i.e., 
u(t, x) - g(v, x). 
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3 Backward diffusion-like PDEs system 



Let us consider again a real-valued function / (f , x) on X R'', with f(t, 0) = 0, which 
has continuous partial derivatives of the first order with respect to f, a - l,m and of 
the second order in x". If 

w, = (w^), t = (fi, r) € Rf, fl = M 

is a multitime vectorial Wiener process, then the function / defines a stochastic process 

3;, = /(f,W,), f eR^ 
By It6-Udri§te Lemma S, llT3l . llT4l . we can write 

dy, = c„ (f) A,/(f, W,) + ^ (t, W,)] df + ^ (t, Wr) dW^, (3.1) 



where 



/32 f 



In order that the stochastic process be a martingale, the drift coefficients in for- 
mula (1) must vanish, i.e., / is a solution of the backward diffusion-like system 

(f, Wt) + \c, it) AJ(t, W>) = 0, for a^h^. (3.2) 
or 2 

Theorem 3.1 The solution f of the backward diffusion-like system (4) depends on the 
point (f ' , f") only through the product of components v = f ' ■ • ■ f™, i.e., it is a function 
of the volume v = f' • ■ ■ f" o/Qq, c R™. 

Proof. From (4) it follows 

cp{t) ^ (f, W,) = Ca(t) ^ (t, W,),a^ p. 

The general solution of this PDEs system is /(f) = ipif ■ ■ ■ f'", W,), f 6 R™. ■ 

Let us take the shape of volumetric features into account, recalling that a volumetric 
function is invariant under the subgroup of central equi-affine (i.e., volume-preserving 
with no translation) transformations, where the determinant of the representing matrix 
is 1. 

If V = f ' ■ • ■ f"' is the volume of Q.q, c R"' and g (v, x) ''^ ip{t^ ■ ■ ■ t"\ x) has contin- 
uous partial derivatives of the first order in v and of the second order in x, then from 
(13.2b it follows that g satisfies the backward multitime heat PDE 

%-\^8-0. (3.3) 
Consequently, the relation (13. Il l reduces to 



dg{v,W,)^^^{v,W,)dW!. 
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Thus, if E [ll^ (v, Wf) Ipj is bounded with respect to f, in bounded subsets of R™, 
then the stochastic process g (v, W,) is differentiable with respect to W" and the par- 
tial stochastic derivative is ^ (v, W,). 

Let T be the maturity multitime of the option. 

Lemma 3.2 Let 

= (W^y, t = (f', f'") eWl, xeR^a^hd 

be a family of Brownian sheets with P{Wf = x) = 1 for t e dR"'. If f : R'^ —> R is a 
continuous bounded function and w(t, x) — EixlfiWj)] is the backward mean value, 
then w(t, x) is a unique bounded solution of backward heat type PDEs 

dw 1 

^(f, J^) = -2 Ca(t)AMt,x), t eRi; w(T,x) = f(x) 
(terminal value problem). 

Proof. Let V - ■ ■ - T'" and v = f' ■ • ■ f™. The transition density is now 
ru ^ 1 l -\\x-y\\^ \ 



Remark 3.3 The heat PDE has two variants, forward and backward, with the forward 
one being more common. 



4 Path independent stochastic curvilinear integral 

The foregoing theory and the papers 0, ||9l, |T3l, fT?! suggest to introduce the notion 
of multitime differentiable stochastic processes. For that we need a multitime vectorial 
Wiener process 

W, = (W^), t = (f', f") € Rf , fl = ITS. 

Definition 4.1 Let jq, c R™ be an increasing curve joining the points 0,t & R"^. A 
multitime stochastic process <l)(f) = <l)(f, Wi), t e W^: is called differentiable with re- 
spect to Wt, on R™, if there exists a multitime vectorial adapted mesurable process 
cf>(t) - ((/)a(t, W,)), t e R'" such that E^\\(p(t)\f-^ is bounded for t in compact sets q/R™ 
and 



0(0 = cD(0) + r 4>a{s)dW:, (4.1) 
•J rot 



as a path independent stochastic curvilinear integral. 

In terms of stochastic differentials, the multitime stochastic process O, is differen- 
tiable if the stochastic system 

JO(0 = Mt) dWf. 
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is completely integrable, i.e., <l)(f) is a function of v = t ■ ■ -f" and hence 0(f) is a 
vectorial function of v [f), ifTSl . 

The multitime process (pa(t) is called the stochastic partial derivative of the multi- 
time process <l)(f) with respect to W" (see also lU, |[3l). 

Remark 4.2 T/ze differentiable multitime stochastic processes have properties simi- 
lar to holomorphic functions: a differentiable process has a differentiable derivative, 
the curvilinear integral primitive of a differentiable process is differentiable, and each 
differentiable process admits a power series expansion. 

Theorem 4.3 Let w(t, x) — E,^y[f(Wj)] be the backward mean value. For any Wiener 
process Y,, the stochastic process w{v, Y,) is differentiable. 

Proof. Based on the previous Lemma, the backward mean value function w(v, Y,) is a 
solution of the backward multitime heat PDE (terminal value problem) 

^ + \Ag^Q, g{T,x)^f{x). (4.2) 
av 2 



5 Hermite polynomials and stochastic processes 

There is a special class of solutions of the backward multitime heat equation which will 
be particularly interesting in stochastic problems. These are the Hermite polynomials 
(see also ifBl ). 

To develope further our theory, we need a point t - (f\ f) € R™, v = f' ■ ■ ■ f", 
a point X - (x', ...,x^) e , a point n - («i, ...,«(/) e N'^, called multiindex, and the 
notations \n\ - n\ + ... + and n! = ni! ■ ■ -ndl. Also for^ = (^', ...,^"') e we write 
^ To simpHfy, the partial derivative 

will be denoted by d^nf. Now, let us use the generating function expansion 

g<x,f>-i,.||fip ^ H„(v,x)^", x,^eR\ ve R^. (5.1) 

neN^ 

Since 

H„(v,x)=l5e..^-'-'f^-^""*|f.o, 
«! 

it follows that H„(v, x) is a polynomial of degree n in jc = (x^, x'^); it is called the n-th 
Hermite polynomial. On the other hand, the generating function g^-'-^^^s' ll^ii" is a tensor 
product of d one-dimensional generating functions of the form e'^ T-^^ . Consequently 

H„{v- x\ x'') = (v; x') ■ ■ ■ H„/v; 
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where H„^ denotes the nk-\h one-dimensional Hermite polynomial. 

Completing the square in the definition of generating function, the n''' Hermite 
polynomial of the variables v e R+ and x e R'' can be written as 

(-V)" Wi _l||,_,,f||2 (-V)" Mi _wd 

H„(v,x) = — —e 2. 5f«e ^" |f=o = — ;— e d^m . (5.2) 
n! n! 

Here and in the rest of the paper, we write v" = v"' ■ ■ ■ y"''. 

The sequence {//„(v, OljtLo ^ complete orthogonal system with respect to the 

weight (2;7rv)^ e'^ . The orthogonality means 

10, \im + n, 
v" (5.3) 
— , it m — n. 
m 

We introduce the multiindex kg - (0, k, 0), with k in the place a. Then, using 
the generating function expansion, taking the partial derivatives with respect to v and 
x", and equating the coefficients, we find 



d 1 d 



dv 2^ dx' 

a 

Consequently each Hermite polynomial i/„(v, x) is a solution of the backward mul- 
titime heat PDE. Thus 

Theorem 5.1 Each Hermite polynomial Hn(y, Wt) is a dijferentiable process and its 
partial derivative with respect to W" is H„-i^(v, Wt). 

It follows that finite sums of Hermite polynomials processes are difFerentiable pro- 
cesses. We extend this statement to series of Hermite polynomials, following the ideas 
of Cairoh and Walsh 0. 

Theorem 5.2 Suppose {a„} is a multiindex family of real numbers such that 
^ fl^ — < oo, for all V > (convergent series). 

Then, the process <!>, defined by 

is dijferentiable relative to Wt and its partial derivative (p,a with respect to W" is given 
by 

4,,„^ ^ a„Hn-iSv,W,). (5.5) 
The convergence of series is understood in L^. 
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Proof. By the orthogonality relations (15.31 1. we find 

^ a„H„ (v, W,) 



v|"l=o 



|«|=0 



This mean value is bounded due to the convergence of the series in the right hand 
member. It follows that the series ( 15.41 ) converges in L? and the same is true for the 
series i5.5i . Consider now the sequence of partial sums 



2 a„H„-iJv,W,), 



|«-1J=0 



and let 



where 



^ \p\ 
•^yot i/i-ui=o 



jI : f' = c\..., f-^ = c"-\f = t" e [0, f] , f^^ = c"^\..., r = c"\ 
is an increasing curve joining the points and t, with = const > Q, [3 + a. It follows 

hm = O, in 

|p|-*oo 

To finish the proof, we need only check that 
lim 



lim (l)^l'^dW"„a^ l,m. 

\p\->co 



Again, by the orthogonaUty relation ( 15.3b . we have 



and consequently 



\n\=\p\+\ 



y«l . . . y«<.-l . . . y"rf 



|n|=|p|+l 



In other words. 



lim r A^i'^/c^ r 
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Theorem 5.3 Suppose that ip (v, x) has continuous partial derivatives of the first order 
in V and of the second order in x, and that [ip (v, Wt)} is a differentiable process, where 
V = f ' ■ ■ ■ f". Then, for each t e R™, we can write 



|«|=0 

2 ^ ^ 1- * — /" f 1 



where the convergence is understood in L and, for t — (r, e R'" 



fl„ = - E[^(v,W,)i/„(y,W,)]. (5.6) 

V 



Proof. Define the process = li|^|=o i^, Wt). Due to the convergence 



< oo, 

" n! 

|«|=0 



the series defining O, converges in and the process O, is differentiable. 
Since //„ is an orthogonal sequence, we get 



; [<D,i/„ (v, W,)] = E [fl„//2 (v, IV,)] = fl„ ^. 



V" 

But, by hypothesis, 



E[^(v,W,)H„(v,W,)] = ^ fl„. 

«! 

Thus, 

E[(p(v,W,) - (t>r) H„(v,W,)] ^ 0, 
for all n € N^. Then, (p (v, W,) - O, = 0. ■ 



6 Union of Tzitzeica hypersurfaces 

A hypersurface M c R™, m > 3, is called Tzitzeica hypersurface Q, lH, provided 
there exists a constant a € R such that we have K - ad'"^^, for all points t = 
(f ' , f") e M, where is the Gauss curvature of the hypersurface and d is the distance 
from the origin of the space to the tangent hyperplane to the hypersurface at the current 
point t. Since the Gauss curvature K describes the shape of the hypersurface, a Tzitze- 
ica hypersurface has a bending against the tangent hyperplane in fixed proportion to the 
normal component of the position vector f. The most simple Tzitzeica hypersurfaces 
are the constant level sets : t^ ■ ■ • t"^ - c in R™ (2™"' connected components). 

Remark 6.1 The Gauss curvature of a Cartesian implicit surface 

Mc : Fit\t^,P) = c 

in M? is 

K = [[F3(F3Fn - 2FiFi3) + f\F^^W3{F^F 22 - 2F2F23) + ^^^33 
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- (F3(-FiF23 + F3F12 - F,Fi3) + FiF2F33f 1 [FjiFj + fI + FlfyK 

The surface Mc is curving like a paraboloid if K(t) > 0, hyperboloid if K{t) < 0, or a 
cylinder or plane ifK(t) = 0, near a point f = t^, E Mc. 

Let us show that the constant level sets of the functions (p(t^ ■ ■■t'", x), with respect 
to t = {t^,...,f) € R™, are Tzitzeica hypersurfaces or unions of simple Tzitzeica hy- 
persurfaces in R™ indexed by the points x. 

Theorem 6.2 If (p : R — >R is a nonconstant function and c is a noncritical value 
of ip, then the constant level set Nc : <pit^ ■■■f) = c is a union of simple Tzitzeica 
hypersurfaces. 

Proof. . Let Ac the set of the solutions of the equation ^(v) = c and v = t^ ■■ ■t"'. Then 
the constant level set Nc is the union of the constant level sets t^ ■■■f" = k, k e Ac. If 
c is not a critical value of <p, then the set Nc is a union of hypersurfaces. ■ 

Remark 6.3 If c is a critical value of ip, then the set Nc is a union of constant level 
sets t^ ■ ■ ■t'" - k, k € Ac, but those level sets which correspond to <p'{k) = are not 
hypersurfaces. 

Theorem 6.4 Each section x = c of the hypersurface H„{y, x) = in W"'^'' = {(t, x)] is 
a union of Tzitzeica hypersurfaces in R+. 

The constant level set 

E[H„(v,W,)H„iv,W,)]=c>Q 
is a simple Tzitzeica hypersurface. 
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